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Color superconductors in which quarks of the same flavor form Cooper pairs are investigated. 
These Cooper pairs carry total spin one. A systematic group-theoretical classiflcation of possible 
phases in a spin-one color superconductor is presented, revealing parallels and differences to the 
theory of superfluid ^He. General expressions for the gap parameter, the critical temperature, and 
the pressure are derived and evaluated for several spin-one phases, with special emphasis on the 
angular structure of the gap equation. It is shown that the (transverse) color-spin-locked phase is 
expected to be the ground state. 
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C ■ It is well-known that certain metals and alloys exhibit a superconducting phase below a critical temperature T^. In 
. this phase, electrons in the vicinity of the Fermi surface form Cooper pairs which leads to a gap in the quasiparticle 
' excitation spectrum 0| . An attractive interaction between the electrons is provided by the exchange of virtual phonons, 
\^ , and the electromagnetic ^^S^ group U{l)em is spontaneously broken. A similar mechanism occurs in sufficiently cold 
' and dense quark matter In this case, the attractive color-antitriplet channel is responsible for the formation of 
\ quark Cooper pairs. Due to the intrinsic properties of quarks (color, flavor, electric charge), many pairing patterns 
' seem to be theoretically possible. In other words, besides the electromagnetic gauge group, also the color gauge group 
C^l \ SU{3)c, the flavor group SU{Nf)f, and the baryon number conservation group U{1)b may be broken completely or 
to a certain residual subgroup. In recent years there have been interesting works studying the ground state of cold 
and dense quark matter, i.e., it has been investigated which color-superconducting phase is favored for certain ranges 
of the quark chemical potential ^. This question is also of phcnomenological interest, since matter in the interior 
of neutron stars can reach densities up to an order of magnitude larger than the nuclear ground state density while 
the temperature can be of the order of keV. Therefore, the core of a neutron star can be expected to be a color 
Q ■ superconductor and the question arises of which of the color-superconducting phases it consists. 
^ ' For very large densities, where the quark masses of the u, d, and s quarks can be considered degenerate since 
• they are much smaller than the chemical potential, the ground state is the so-called color-flavor-locked (CFL) phase 
0. In this phase, quarks of all flavors and all colors form Cooper pairs, breaking spontaneously SU{3)c x SU{3)f x 
U{l)em X U{1)b to a subgroup SU{3)c+f x U{l)c+em- At smaller densities the situation is more complicated, because 
' the strange mass cannot be neglected. Furthermore, the conditions of /3-cquilibrium and electric charge neutrality 
^ [ impose restrictive conditions on the system. The simplest solution seems to be a pairing of only u and d quarks Q], 
the so-called 2SC phase. In this phase, quarks of one color remain unpaired, and the symmetry breaking pattern is 
SU{3)c X SU{2)f X U{l)em X U {!) B ^ 5'C/(2)c X SUi2)f x U{l)c+em X Uil)em+B- In both CFL and 2SC phases the 
Fermi momenta of the quark flavors participating in pairing are assumed to be equal. This is a necessary condition for 
the conventional BCS pairing mechanism, since in both cases quarks of different flavors form Cooper pairs. However, 
for moderate densities, this assumption is not valid. Therefore, the ground state is neither the (pure) CFL nor the 
(pure) 2SC state. Several other possibilities have been discussed. In principle, besides a transition to the normal- 
conducting state, there are two classes of alternatives, both yielding color-superconducting states. In the first class, 
the difference in Fermi momenta is nonzero but small enough to still allow for pairing of quarks of different flavors. 
The second class accounts for cases in which quarks of the same flavor pair. 

Let us first mention some options for the first class. First, there might be a phase in which the Cooper pairs 
carry nonzero total momentum. In this case, the system exhibits a crystalline structure due to a spatially varying 
energy gap Q . This kind of superconductivity is called LOFF phase and was first discussed in solid-state physics 1^ , 
where a difference in the electron Fermi momenta is induced by an external magnetic field. Second, there are studies 
about the "gapless" 2SC [H and CFL |3| phases. In these phases, at least a part of the quasiparticle excitations are 
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ungapped, although the gap parameter is nonzero. This feature obviously can have enormous physical consequences, 
for instance for the specific heat and the neutrino emissivity which both affect the cooling of a neutron star with a 
core in a gapless color-superconducting phase. Other possible phases are derived from the CFL phase and contain 
kaon and/or eta condensates Moreover, besides a displacement (LOFF phase), also a deformation of the Fermi 
surface has been discussed 0. 

In this paper, we discuss the possibility of the second class, namely systems in which quarks of the same flavor 
form Cooper pairs 0, llCt llll IT^ . The simplest situation for this kind of color superconductivity is a system of only 
one quark flavor. More realistic scenarios are two- or three-flavor systems where each quark flavor separately forms 
Cooper pairs. Another possibility is a system where u and d quarks form a (gapless) 2SC superconductor while the s 
quarks pair separately. 

Since the attractive interaction of quarks is provided in the anteymmetric color antitriplet channel, the spin channel 
must be symmetric in order to ensure the overall antisymmetry of the Cooper pair wave function. Consequently, 
Cooper pairs consisting of quarks of the same flavor cannot carry total spin zero but must condense in the spin-one 
channel, J = 1, where J ~ L + S is the total spin of the Cooper pair, consisting of their angular momentum L and 
their spin 5*. This "two-triplct-condensation" (color and spin triplets) renders the structure of the order parameter a 
complex 3x3 matrix. This is in contrast to the 2SC case, where, due to condensation in the color antitriplet, flavor 
and spin singlet channels, the order parameter is a complex 3-vector. In the CFL case, the order parameter is also a 
3x3 matrix, originating from the color and flavor antitriplets. 

Another system with this structure of the order parameter is superfluid ''He In this nonrclativistic system, 

angular momentum L and spin S are separate quantum numbers, both giving rise to a triplet structure of the 
condensate. In other words, a Cooper pair of '^He atoms carries angular momentum one and spin one. Without 
external magnetic fields, two phases of superfluid ^He are experimentally known, the A and B phases. The A phase 
is given by the order parameter structure Ay = 5i'i{5ji + iSj2), where the index i refers to spin and the index j to 
angular momentum. This order parameter spontaneously breaks SO{3)s x SO{3)l x U{1)n to U{l)s x U{1)l+n, 
where SO{3)s and SO{3)l are the spin and angular momentum groups, respectively, and U{1)n is the particle number 
conservation group. The B phase, which covers the largest region of the phase diagram, is given by Aij = Sij, locking 
angular momentum with spin. Here, the residual group is SO{3)s+l- While the gap function (^(k) is isotropic in the 
B phase, it is anisotropic in the A phase. The condensation energy (density) Ap of the superfluid states is given by 
the angular average of the square of the gap , 



Ap=-iV(O)(02(k))t, (1) 



where A^(0) is the density of states at the Fermi surface and {—)^ = J ^f^- In weak coupling and at zero temperature, 
the ratio of the condensation energies of the A phase and the B phase is (cf. Eqs. (3.72) and (3.75) of Ref. [Tjl) 

ApA 



ApE 



0.88 . (2) 



Therefore, weak-coupling theory predicts the B phase to be the favored state. 

It is the main goal of this paper to determine the favored state in a spin-one color superconductor. The paper is 
organized as follows: In Sec. ^ we discuss possible symmetry breaking patterns in a spin-one color superconductor. 
In general, the group G = SU{3)c x SU{2)j x ?7(l)em x U{1)b is spontaneously broken to a residual group iJ, where 
SU{2)j is the (relativistic) spin group. A systematic list of order parameter matrices A and the corresponding residual 
groups is presented, based on the simple group-theoretical condition that A be invariant under transformations of 
H. Four phases with large residual groups are picked for further investigation in the next sections, namely the polar, 
planar. A, and CSL (color-spin-locked) phases. 

After establishing the formalism in Sec. IIII Al Sec. IIII Bl is devoted to solving the gap equation for the general 
case of one or two nonzero (constant or angular-dependent) energy gaps. The results are the gap functions for zero 
temperature, T = 0, at the Fermi surface 

0r(k) = \/%^(f>o , (3) 

where the angular dependence is contained in the quantities Afc^^ = Ar(k). The gap function occurs in the excitation 
spectrum of the quasiparticles. 



y^(fc-/i)2-|-Afc.,02 (4) 
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Let us briefly recall the situation of isot ropic gap f unctions, Xk._r = Ap. It has been shown that, in this case and with 
r — 1,2, the gap parameter is given by 



c ropic gap t unctions, A 

111 in mill III 111 



o = 26&;,e-''e-Vexp(-^) , (5) 



2^; ' 



where 



' ^ '''''' iJ^T ^ 

with the strong coupling constant g. The exponent d is zero in all spin-zero phases and nonzero for J = 1. The 
exponent ^ is defined as 

C^ln(ArA-)l/^ (7) 

where the numbers ai and 02 have to be determined for each phase separately. They fulfil the condition 

ai + aa = 1 . (8) 

The notation e^^ instead of (A°^ X'^^)^^/'^ will turn out to be convenient for the generalization to the case of anisotropic 
gap functions. 

In the case of the 2SC phase, there is only one gapped quasiparticle excitation branch, ai = Ai = 1, 02 = A2 = 0, 
and hence 

cf>f^ = 2hh'^^xe^^{-^^ , (9) 
whereas in the CFL phase, there are two different gaps, ai — 1/3, 02 = 2/3, Ai = 4. A2 = 1, which leads to 

^CFL ^2-1/3 02SC^ 

In Sec. nil CI we compute the transition temperature Tc for the transition from the normal-conducting to the 
superconducting state. As for the gap, this is a generalization of the cases with constant (= angular-independent) 
gaps, for which the transition temperature is given by (19| 

Z£ = £le^~0.57e« , (10) 

where 7 ~ 0.577 is the Euler-Mascheroni constant. It has been one of the main conclusions of Ref. il9|| that this 
expression shows the violation of the well-known BCS relation Tc/(t>o — 0.57 in the case of a two-gap structure, i.e., 
Ai,2, ai,2 7^ 0. (Also in the case of the gapless 2SC phase, this relation is violated 0.) 

As in the case of '^He, one expects the preferred phase to have the largest condensation energy, cf. Eq. (Q. We 
specify this statement in Sec. IIII Dl with a general derivation of the pressure in an arbitrary color-superconducting 
phase. 

In Sec. IIVI we determine the excitation spectrum, the gap functions, the critical temperature, and the pressure 
for the polar, planar. A, and CSL phases. For each phase we consider three special cases, termed "longitudinal", 
"mixed", and "transverse". These three cases arise from the following property of spin-one phases: Contrary to a 
spin-0 color superconductor, where only quarks of the same chirality form Cooper pairs (RR and LL pairs), in a 
spin-1 color superconductor also pairing of quarks with opposite chirality (RL and LR pairs) is possible |lll ll2L 
In general, the order parameter contains a linear combination of both kinds of condensates. We focus on the cases of 
pure RR/LL condensates (longitudinal), a special admixture of RR/LL and RL/LR condensates (mixed), and pure 
RL/LR condensates (transverse). Consequently, in total 12 phases are studied. 

In Sec. we summarize the results and give an outlook for possible consequences in neutron stars with color- 
superconducting cores. 

Our convention for the metric tensor is g^'^ = diag{l, —1, —1, —1}. Our units are h = c = ks = Four- vectors are 
denoted by capital letters, K = K'^ = {ko, k), and k = |k|, while k = k/fc. We work in the imaginary-time formalism, 
i.e., T/VJ2k = ^X^n / '^'^k/(27r)'^, where n labels the Matsubara frequencies w„ = iko. For bosons, ujn = 2miT, for 
fermions, w„ — (2n + l)7rT. 
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II. PATTERNS OF SYMMETRY BREAKING 



In this section, we discuss possible symmetry breaking patterns in a spin-one color superconductor. In other words, 
we present a systematic classification of theoretically possible superconducting phases. In the case of a one-flavor, 
spin-one color superconductor, the relevant original symmetry group of the system is 

G = Gi X G2 X G3 , (11) 

where 

Gi = 5(7(3), , G2 = SU{2)j , G3 = UHU, . (12) 

Note that the global group U{1)b, accounting for baryon number conservation, has the same generator as the local 
symmetry group U{l)em- This is not true in a system with Nf > 1, when at least two quark flavors differ in 
their electric charge. In particular, in the CFL phase {Nf — 3), this leads to the fact that the system is not an 
electromagnetic superconductor but a superfluid (the breakdown of U{1)b gives rise to a Goldstone boson). For 
Nf = 1, however, these two phenomena are coupled, i.e., a superflow is equivalent to a supercurrent. 

The order parameter A is an element of a representation of G. In the following, we use the term order parameter 
somewhat sloppily for the pure matrix structure A. For a spin-one color superconductor, the relevant representation 
of G is the tensor product of the antisymmetric color antitriplet [3]JJ and the symmetric spin triplet [3] j, 

A e [3]? ® [3]} . (13) 

Therefore, A is, as in the case of superfluid '^Hc, a complex 3x3 matrix. There is no nontrivial contribution from the 
flavor structure since we consider systems with only one quark flavor. The group G is spontaneously broken down to a 
residual (proper) subgroup H Q G. This means that any transformation g (z H leaves the order parameter invariant, 

.9(A) = A . (14) 

In the following, we investigate this invariance condition in order to determine all possible order parameters A and 
the corresponding residual groups H. The method we use in this section is motivated by the analogous one for the 
case of superfluid '^He 14] . 

First, one has to specify how G acts on the order parameter in Eq. (|14|l . To this end, we write an arbitrary group 
element g G G as (in this section, no confusion with the strong coupling constant g is possible) 

9 = (51,52,53) , (15) 

where 

gi = exp(-ia,„r^) , 52 = exp(j5„J„) , 53 = exp(2jcl) , (16) 

with real coefficients a„i {m = 1,...,8), 6„ {n — 1,2,3), and c. The Gell-Mann matrices Tm generate the group 
5t/(3)c, and we have taken into account that the color representation is an anMriplet. The matrices J„ are the 
generators of the spin group SU{2)j, {Jn)ij = —i f-nij- For the generator of U (1) we choose 2 • 1, where 1 is the 3x3 
unit matrix. The factor 2 accounts for the diquark nature of the order parameter. 

Let us now introduce a basis Ji ^ kj {i,j = l, 2, 3) for the representation given in Eq. (|13|l . For the basis elements 
of [3]° we choose the antisymmetric 3x3 matrices Ji, as introduced above as generators of the spin group. The basis 
of [3]j is given by the 3-vector Kj, which will be specified in Sec. IIVI cf. Eq. H12U|) . Thus, we have to consider the 
action of G on 

TWk ^ J^ Ay-Kj . (17) 

We have 

g{M^) = g3g{''JkA,,gi'^i. (18) 

Therefore, the matrix A transforms as 

5(A,,)=g3 5f A«5^^ . (19) 
Then, using Eqs. H16|) . the infinitesimal transformations of A by G are given by 

5(A) ~ A - a^T^A + 6„ A J„ + 2cA , (20) 
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where T^A as well as AJ„ are matrix products. The invariance condition for the order parameter H14() is thus 
equivalent to 

-amT^A + 6„AJ„ + 2cA = . (21) 

This matrix equation can be written as a system of nine equations for the nine complex entries An, ... , A33 of the 
matrix A. In principle, one can find all possible symmetry breaking patterns and corresponding order parameters by 
setting the determinant of the coefficient matrix to zero. Then, each possibility to render the determinant zero yields 
a set of conditions for the coefficients a^, 5„, c, and it can be checked if these conditions correspond to a residual 
subgroub H . But since this is much too complicated, we proceed via investigating possible subgroups explicitly. In 
the following, we focus only on the continuous subgroups of G. 

Let us start with subgroups H that contain the smallest possible continuous group, U{1), i.e., 

H = [/(I) X H' , (22) 

where H' is a direct product of Lie groups. The residual U{1) must be generated by a 3 x 3 matrix U which is a linear 
combination of the generators of G, i.e., in general, 

U = amTm + hnJn + 2cl . (23) 

Let us restrict ourselves to linear combinations that involve one generator of each group Gi , G2 , G3 , for instance 

C/ = a8T8 + fe3J3 + 2cl . (24) 

With Eq. (|24|l . the invariance condition 

e*^(A) = A (25) 

results in a system of nine equations, which can be discussed explicitly. This is done in Appendix We find 
ten different order parameters, eight of them depending on two or more parameters Ai,A2,..., which physically 
corresponds to more than one gap function. With the normalization 

Tr(AA^) = 1 , (26) 



the calculation in Appendix^yields also two matrices A that do not depend on any free parameter (which corresponds 
to only one gap function). In analogy to "^He, let us call the phases defined by this kind of order parameter inert 
All experimentally known states of superfluid ^He belong to this class of order parameters. Mathematically speaking, 
these matrices play a special role due to a theorem ( "Michel's Theorem" ) 0, , which ensures that these order 
parameters correspond to a stationary point of any G-invariant functional of A (for instance the effective potential). 
This is the reason why we also focus on these order parameters in the explicit calculations of the physical properties, 
see Secim 

One order parameter, found with the ansatz H24|l and corresponding to an inert phase, is 

(27) 

defining the polar phase. It has its analogue in '^He [T^ . where the order parameter matrix is identical. For the case 
of a color superconductor, certain aspects of the polar phase have already been discussed in Refs. 0,0,0,112. For 
the corresponding residual group H see Fig.^ where we list all order parameters. More details, especially the explicit 
forms of the generators of H, are given in Appendix^ The second order parameter giving rise to an inert phase is 








1 i 



A = I I , (28) 



This order parameter leads to the A phase |0, Q| . 

In order to find all possible (inert) order parameters, it is necessary to consider at least one more combination for 
the residual U{1) in Eq. (j^ . namely 



U ^a2T2+b3J3 + 2cl 



(29) 
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The reason why T2 plays a special role is that we used the generator J3 of the spin group, which is proportional to T2. 
Consequently, we expect to find additional residual groups that connect the color group with the spin group (meaning 
a residual U{\) generated by a combination of a color and a spin generator). The calculations with this generator are 
completely analogous to the ones with the ansatz (|24|l . Therefore, we present the result without elaborating on the 
details. The ansatz (|29|) yields one inert order parameter that is different from the above ones, namely 




(30) 

This order parameter corresponds to the planar phase The residual group is H = U{1) xU{l), with generators 



c/ = 2T2 + J3 , v^n + ^i. (31) 

Let us now tiu'n to possible groups H that do not contain any U{1) but solely consist of higher-dimensional Lie 
groups, say 

H = SU{2) X H' . (32) 

Let Ui, U2, U3 be the generators of the residual SU{2). They are linear combinations of the generators of G, 

= <,r,„ + 6; J„ + 2c'l , i = 1, 2, 3 . (33) 

Since they must fulfil the SU(2) commutation relations, 

[Ji,Jj]=ieijkJk, i,j,k<3. (34) 

they must not contain the generator of G3 = U{1), the unit matrix, i.e., = = c"^ = . Therefore, there are three 
possibilities. First, each Ui is a combination of color and spin generators. Second and third, each Ui is composed 
solely of color or spin generators, respectively. The simplest options to realize these cases are 

h , (35a) 

(35b) 
(35c) 

where (Tl,T^,T^) is either given by (Ti,r2,T3) or (2r7, -2r5, 2T2), which both fulfil the required commutation 
relations. Using the options H35a|l - (|35c|l . let us first show that H' in Eq. (|32|) cannot be a second SU(2). To this 
end, assume that H' = SU{2) with generators Vi, V2, V3, which have the same form as the generators Ui in Eqs. H35|l . 
Then, since the Lie algebra of _ff is a direct sum of the constituent Lie algebras, we have to require 

[U,,V,]^0. (36) 

This condition reduces all options to one, namely 

U^^Tl , V, = J, (37) 

(or vice versa). However, now the invariance equation for the order parameter yields 

A J, = , (38) 

for all 1 = 1,2, 3, which does not allow for a nonzero order parameter A. Therefore, H' = SU{2) is forbidden. Since 
the cases with H' = U{1) and H' = U{1) x U{1) were already covered in the above discussion, the only possibility 
that is left is = 1 and thus H = SU{2). 

The generators in Eq. (|35c|) can immediately be excluded since they also lead to Eq. (|38|l . The same argument 
excludes case pSKj) with (r{,T^,T^) = (2r7, -2T5, 2T2). Case (|35H|) with {T{,T^,TQ = {Ti,T2,T3) leads to two order 
parameters already considered above, namely the polar phase, Eq. H27|l . and the A phase, Eq. (|28|l . In case (|35a|l . 
only {Ti,T^,TQ = (2r7, -2T5, 2T2) is possible. With 





= T 




= T' 







U,^ = -T/A + A J, = , 



(39) 
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one finds 

1 / 1 0\ 

(40) 

Indeed, it can be cliecked witli Eq. (|21|l tliat tliis order parameter leads to 

ai = as = tti = ttQ = as = c = , 02 = 263 , a^ ^ -2^2 , 07 = 2bi , (41) 

wlrich corresponds to H — SU{2), consisting of joint rotations in color and spin space. This is the color- spin-locked 
(CSL) phase, discussed for a spin-one color superconductor in Refs. [1, [la Ull l23l ■ I* is the analogue of the B 
phase in superfluid '^He. 

Finally, we give an argument why an even larger subgroup, i.e., an SU{i), cannot occur in the residual group H. 
Assume that there are eight generators Wi, . . . , of this SU{?>). Then, as for the SU{2) subgroup above, there can 
be no contribution to Wi, . . . , Ws from the G3 generator due to the SU{3) commutation relations for the generators, 

m,W,] = I f,jkWk , *,j,fc<8, (42) 

where fijk are the SU (3) structure constants. Also, Wi — Ti is excluded because in this case the invariance condition 
yields TgA = and thus A = 0. Therefore, at least one of the spin generators has to be included. For instance, 
choose Ws = Tg + 63 J3. Then, from the commutation relation [W4, W5] ~ W^, we conclude that also Ji and J2 must 
be included via W4 = 14-1- 61 Ji and = T5 + &2>/2. But now the three equations W4A = W5A = Wg^ = lead to 
A = 0. Therefore, we conclude that there is no residual group H that contains an 5[/(3). For a more rigorous proof 
one has to take into account more complicated linear combinations of the original generators. 

In Fig.^ we summarize our results in a list of all superconducting phases that we have found in the above discussion. 
It should be mentioned that this list is not complete, since for the generators of the residual C/(l)'s we have restricted 
ourselves to two special forms given in Eqs. (|24|l and H29() . Therefore, there are certainly more (at least non-inert) 
order parameters that lead to an allowed symmetry breaking. The inert states are listed in the first four lines. Each 
of these four states has its analogue in superfluid ^He. Note that the Ai phase, which is experimentally observed in 
■^Hc in the presence of an external magnetic field, does not lead to an allowed symmetry breaking in the case of a 
spin-one color superconductor p^. To see this, one inserts the order parameter of the Ai phase, 

(43) 




into Eq. H21|l . One obtains 

ai = as = ^1 = ^2 = , a4 = ay , 05 = ag , -02 H ^ag — 63 — 2c = . (44) 

^ 2v 3 

These 7 conditions suggest that dimi? = 12 — 7 = 5 and thus H — SU{2) x U{1) x U{1). However, there is no 
possibility to construct three generators from the above conditions that fulfil the SU{2) commutation relations. For 
instance, assume that two of these generators are given by Ui = T4 + and U2 = T^, + T^. Then, with Eq. 142|l . 
[Ui, U2] ~ T3. But since 03 = 0, the third generator U3 cannot be proportional to T3. Consequently, there is no Ai 
phase in a spin-one color superconductor. 

Below the four inert states we list the eight non- inert states which have been found in Appendix ^ and which we 
term Pi, . . . ,Pg. Note that one of these non- inert phases. Pi, has a larger residual symmetry group than the planar 
and CSL phases. 

There are several properties of the spin-one phases which can easily be read off from Fig.^ First, consider the spin 
group SU{2)j. This symmetry accounts for the rotational symmetry in real space of the normal-conducting phase 
(in the case of the spin-one representation, one can equivalently consider S0{3)j instead of SU{2)j). The list shows 
that spatial symmetry is broken in each case. For instance, in the polar phase, SU(2)j is broken to its subgroup 
U{l)j. Therefore, the superconducting phase is invariant under rotations around one fixed axis in real space. In most 
of the other cases, the breaking of the spatial rotation symmetry is more subtle: For instance in the planar phase, 
the superconducting state is invariant under a special joint rotation in color and real space. The most interesting 
breakdown of spatial symmetries is present in the CSL phase. Here, any rotation in real space leaves the system 
invariant as long as one simultaneously performes the same rotation in fundamental color space which is spanned by 
the three directions red, green, and blue. 
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FIG. 1: Symmetry breaking patterns and order parameters for a spin-one color superconductor. The original symmetry group 
(first line) is given by the color gauge group SU{3)c (blank background), the spin group SU{2)j (grey background), and the 
electromagnetic gauge group U{l)em (hatched background). The backgrounds of the residual groups illustrate the symmetry 
breaking pattern. For instance, the blank SU{2) occuring in the residual group of the polar phase is generated solely by 
generators of the original SU{3)c while the blank/hatched f/(l) in the same line is generated by a linear combination of the 
generators of the original SU{2)j and U{l)em groups. For the explicit expression of these generators, see text. 



Next, let us read off some properties concerning the color symmetry. It is obvious that in none of the cases the 
full color symmetry is preserved. In this sense, it is justified to call each phase a color superconductor. In three 
of the cases, there is a residual color subgroup SU{2), namely in the polar phase, the A phase, and the Pi phase. 
Mathematically speaking, this residual group originates from the fact that the order parameter has only nonzero 
elements in its third row. Therefore, the third direction in fundamental color space is preferred. Physically, this 
means that the Cooper pairs carry color charge anti-blue, or, in other words, only red and green quarks form Cooper 
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pairs. Of course, the choice of the anti-blue direction is convention; more generally speaking, quarks of one color 
remain unpaired. Remember that this is also true for the 2SC phase. 

A spontaneously broken gauge symmetry gives rise to massive gauge bosons. In the case of a color superconductor, 
these masses are the magnetic screening masses of the gluons. Therefore, in the cases where there is a residual color 
subgroup SU{2), we expect a Meissner effect for five of the eight gluons. Three of the gluons, however, namely those 
corresponding to the generators Ti, T2, T3, do not attain a Meissner mass. This is also obvious from the fact that 
these gluons do not see the (anti-)blue color charge which is carried by the Cooper pairs. Also with respect to the 
breakdown of the color symmetry, the CSL phase is exceptional. Although there is a residual SU{2), where three of 
the color generators are involved, we expect all eight gluons to attain a Meissner mass. To this end, note that this 
residual SU (2) is a global symmetry and therefore all dimensions of the gauge group have to be considered as broken. 
This is analogous to the CFL phase, which also exhibits a color Meissner effect for all eight gluons [sj. For a more 
detailed and quantitative discussion of the color Meissner effect, see Ref. [2^. 

In order to discuss the question whether the color superconductors in Fig. ^ are also electromagnetic supercon- 
ductors, one has to consider the hatched backgrounds in the residual groups. From ordinary superconductors, we 
know that U{l)em is spontaneously broken below the critical temperature. Therefore, a simple conclusion is that all 
states in the list without a hatched background occurring in the residual group are electromagnetic superconductors. 
Obviously, this is the case for the CSL phase and for one of the non- inert states, namely the Pk phase. The electro- 
magnetic Meissner effect in the CSL and polar phases has been discussed in detail in Refs. |2lL 12^. In particular, 
it has been shown that also the polar phase exhibits an electromagnetic Meissner effect in the case of a many-flavor 
system, where at least two of the quark flavors have nonequal electric charges. 



III. GAPS, TRANSITION TEMPERATURE, AND PRESSURE 



In this section, we derive general expressions for the gaps (Sec. IIIIB|I . the transition temperature (Sec. IIII Cp . 
and the pressure fSec. IIII D|l . valid for all color-superconducting phases we consider. In Sec. II VI we evaluate these 
expressions for several spin-one color superconductors. 



A. Notations and definitions 



The starting point is the effective action F, which yields the gap equation as well as the effective potential VcS = —y^ 
and the pressure p, which is the negative of the effective potential at its stationary point, p = —VcS- 

The effective action F can be derived from the QCD partition function using the Cornwall-Jackiw-Tomboulis (CJT) 
formalism '23^ . The resulting functional can be written as |3j HM US 



1 



1 



1. 



1. 



T[Dg,Df] = --TrlnD^' - -Tv{A-' Da - l) + xTrlni^^^ + -Tt{S^'Df - 1)+T2[Dg,Dp] . 



(45) 



Here, Dq (Aq) and Dp (Sq) arc the full (tree-level) gluon and fermion propagators, respectively. All propagators are 
defined in Nambu-Gor'kov space. We account for the doubling of degrees of freedom in this basis by introducing the 
factor 1/2. The traces run over Nambu-Gor'kov, Dirac, flavor, color, and momentum space, and T2[Dg, Dp] denotes 
the sum of all two-particle irreducible diagrams. In the following, we will consider the two-loop approximation of this 
sum which, for the fermionic degrees of freedom, is equivalent to taking into account only one diagram (the "sunset" 
diagram 26]). The stationarity conditions for T[Dg^ Dp] yield Dyson-Schwinger equations for the inverse gluon and 
fermion propagators. 



= S,-' 



n , 



(46a) 
(46b) 



Here, the pair of propagators (A, 5) is the stationary point of the effective potential and we defined the gluon and 
fermion self-energies as 



n = -2 



ST. 



SDg 

The free inverse fermion propagator is 



E = 2- 



STn 



{Dg,Df) = (A,S) 



5Di 



{Dg,Df) = (A,S) 



So' 



[Go]" 



(47) 



(48) 
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where [Gq] ^ and [Gq ] ^ are the free inverse propagators for massless quarks and massless charge-conjugate quarks 
in the presence of a chemical potential fi, 

[G±]-i=7^^M ±M7o . (49) 

In order to find the full propagators, one has to solve the Dyson-Schwinger equations, Eqs. (|46a|l . self-consistently. 
To this end, we denote the entries of the fermion self-energy in Nambu-Gor'kov space by 

^^(T+ t-) - (50) 



and invert Eq. (|46b|l formally p^ . which yields the full quark propagator in the form 

5=(£^:), (51) 

where the fermion propagators for quasiparticles and charge-conjugate quasiparticles are 

G± = {[G±]-i + I]± - $T([Gj]"i -f , (52) 

and the so-called anomalous propagators, typically nonzero for a superconducting system, are given by 

S± = -([G^]-i + ST)-i$±G± . (53) 
In the two- loop approximation of the quark self-energy S(_ftr) in momentum space is 

m) = E s{Q) n K'uiK ~ Q) , (54) 

Q 

where 

-1' 



'q" ^StT ) ^ (55) 



with the Gell-Mann matrices Ta, a — 1, ... ,8. Due to the Nambu-Gor'kov structure, Eq. (|54|l is actually a set of four 
equations. With Eqs. (|5()|l and (plj, the off-diagonal (21)-component leads to the gap equation 

^+{K)=g'^Y.^^T;[E+iQ)rn A^t (K^Q). (56) 
Q 

The quantities ^^(K) are matrices in flavor, color, and Dirac space and functions of the quark four- momentum K. 
Both quantities are related via 

*"=7o[*+]S. (57) 

Following Ref. |l9l |. we term the matrix 4>+(if) gap matrix and use, for condensation in the even-parity channel and 
in the ultrarelativistic limit, the ansatz 

^+{K) = Y,MK)Mi,Ai. (58) 

e=± 

The Dirac matrices = (1 + 6797 • k)/2, where e — ±, are projectors onto positive and negative energy states, and 
4>e{K) is the gap function. The quantity is a matrix in color, flavor, and Dirac space. It is defined by the order 
parameter and thus determines the color-superconducting phase. For the explicit form of A^k in the case of a spin-one 
color superconductor, see Eq. (|120|l . We can always choose Mk such that it commutes with the energy projectors, 

[A^k,A^]=0. (59) 

The diagonal elements of the quark self-energy can be approximated as ^3 

S+ =S- ~7off'fco In-rr . (60) 

Kq 
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where = (37r/4)TO^; the zero-temperature gluon mass parameter (squared) is = N f ^"^ / [QTr"^) . 
Usmg Eqs. (|52|) . H58|) . and (|6U|) . we find for the fermion propagators 

where 



Z(fco)^ (l+g^n-^) (62) 



is the wave function renormahzation factor introduced in Ref. 27]. In Eq. (|61|) we introduced two sets of projectors, 
V^^ and "Pkr! they project onto the eigenspaces of the matrices 

^^ = 7oXtXk7o and L- = M^mI , (63) 
respectively, i.e., if the number of different eigenvalues \k,r is the projectors are |28| . 



In the case of spin-zero color superconductors (2SC and CFL), and are identical. But this is not true for all 

spin-1 phases as we shall see in Sec. IIVI However, in all phases we consider, and have the same spectrum. 
Consequently, 

Li ^Y.>'^^-Kr ■ (65) 

r 

We denote the degeneracy of the eigenvalue Afc,r by 

^ TriT'^J . (66) 

In general, the eigenvalues Xk,r depend on the direction of the quark momentum k (they do not depend on the modulus 
k). They enter the quasiparticle excitation energies introduced in Eq. (|61|l . 



^{ek-ixY + \k,r\c^e? ■ (67) 



Consequently, the spectrum of the matrices determines the structure of the quasiparticle excitations. In all phases 
we consider, there are at most three different eigenvalues, and at most two different nonzero eigenvalues A^.^. A zero 
eigenvalue corresponds to an ungapped excitation branch. This is for instance the case in the 2SC phase, where the 
ungapped blue quarks give rise to a zero eigenvalue A2 = 0, while Ai = 1. Two different nonzero eigenvalues correspond 
to two excitation branches with different gaps, well-known from the CFL phase, where there is a quasiparticle singlet 
with gap 2(j) (Ai — 4 with degeneracy 1) and a quasiparticle octet with gap (/) (A2 — 1 with degeneracy 8). In the case 
of the spin-1 phases, the eigenvalues carry the potential angular dependence of the energy gap. In the following we 
assume that there is no additional angular dependence in the functions <f>e- Moreover, since we neglect the antiparticle 
gap, 4>- ~ 0, we may denote the particle gap by ^ = 0+. We assume this function to be real, |0p = (j)^, and denote 
the value of this function at the Fermi surface by ipQ. With these assumptions and notations we can define the 
root-mean-square (= quadratic mean) of the function yjXk.r <Po s-s 



'{Xk^rh^bo. (68) 
Furthermore, for the following it is convenient to define the normalized eigenvalue 

A., ^ . (69) 

We make use of these definitions in the final results, cf. Eqs. (|98|) . IjlOSI) and H118() . 

Finally, we determine the anomalous propagator S+ . Inserting the expression for the propagator , Eq. (|61|l , into 
the definition l)53|l and using the form of the gap matrix given in Eq. 1)58(1 . we obtain 

S+(i^) = - E 70 A^k 70 ^k^.Ar , . • (70) 

Fo/^(fco)] "(ew) 
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B. Solution of the gap equation for an anisotropic gap 



In this section, we solve the gap equation for the general case of an anisotropic gap function. Formally, with the 
definitions of the previous section, this means that we allow for angular-dependent eigenvalues A^^^. Starting from 
the gap equation (|56|l for the gap matrix ^^(K), we obtain a gap equation for the function (p^iK) by inserting Eqs. 
(|58|l and H70() into Eq. (|56|l . multiplying both sides with Alj(.A^, and taking the trace over color, flavor, and Dirac 
space. Moreover, we use the gluon propagator in the Hard-Dense-Loop (HDL) approximation, which is permissible to 
subleading order We denote quasiparticle energies by e^ s = e^^, because, due to ~ 0, the quasioniiparticle 
energies do not occur in the gap equation. Since the effect of the wave function renormalization factor Z{ko) for the 
gap is known [l7lll8|| and does not affect the angular structure of the equation, we omit it for simplicity. The resulting 
factor 6q, see Eq. © can easily be reinserted into the final result. We arrive at 



4 EE 



0(Q) 

1o - 4, 



A^''(if-Q)T;,(k,q) , 



(71) 



where, following Ref. [I^, we have defined 



^A(k,q) 



Tr 



7^ T,i' 70 Mq 70 V+, Aq 7. T,,. mI A+ 



Tr 



(72) 



With P = K — Q, the gluon propagator in pure Coulomb gauge is given by 

A0"(P) = A,(P) , AO^(P) = , A'^P) = {S''^ -ffP) At{P) 



(73) 



For the definition of the longitudinal and transverse gluon propagators Ag^t see for instance Ref. With the 

(negative) transverse projection of the tensor "/^"^(k, q). 



(74) 



we obtain after performing the Matsubara sum (for details see Ref. 11]), and after taking the thermodynamic limit, 
V ^00, 



r 
4 



(27) 



3 



f±l tanh ( ^ 



V 2T 



Ft{p) rot,(k, q) + Ft{p, eq.., ek,r) r/(k, q) 



(75) 



where 



arises from static electric gluons, while 

2 



Ft{p, eq,s,ek,r) 



p^ 



Flip) 



Q{p - M) + Q{M - p) 



3m2 



(76) 



(77) 



originates from non-static and almost static magnetic gluons. In Eq. (|75|l we have abbreviated 4>k,r = 0(ek,r,k) 
and 4>q^s = 0(eq.s,q)- At this point, the angular integral (Klq on the right-hand side of Eq. (|75|l seems to be too 
complicated, since, besides the square of the gluon 3-momentum p^ = fc^ + — 2 q • k and the functions T^q, ^ also 
the excitation energies Cq.s depend on the direction of q. The solution to this problem is to multiply both sides of 
the equation with 



4 = Tr 



XkX[A+ 



(78) 



and take the angular average over k of the whole equation. The right-hand side of Eq. I)78|l is obtained with the help 
of Eqs. (|63|l . 1)65(1 . H66() . and the following identities. 



TriP+J = Tr[P+ A£] = Tr[P. ^A^ 



(79) 
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They will become obvious in Sec. IIVI where we discuss the specific phases in detail. The only nontrivial phase with 
respect to these identities is the A phase, where V^^^ ^ ^kr- ^''^ ^^^^ case, one uses Eqs. (|B13|I from Appendix IBl to 
prove these relations. 
We arrive at 

(4)fc0fe.. = ^ 1(03 E^tanh(^) ^F,(p)|To«o(k,q)+^^,b,.,,„.k,.)|7i^(k,q)^_ , (80) 



where we pulled the factor £q out of the k integral on the right-hand side of the equation. This turns out to be 
convenient for the calculation. On the left-hand side we have made use of the assumption that the function (l)k,r does 
not depend on the direction of k. We shall see below that this assumption is consistent with our final result. 

Now the angular integral over k has to be performed for each phase separately. In Appendix |21 we present this 
calculation explicitly for the transverse polar, planar, and A phases, and in Appendix IdI it is presented for arbitrary 
longitudinal phases. However, we can give a general result similar to the method introduced in Ref. 19]. For the 
angular integral we use a frame with z'-axis parallel to q. However, besides the fixed direction given by q, in 
all nontrivial cases there is another preferred direction in the system, namely the one that is picked by the order 
parameter. This has to be taken into account for the integral and is explained in detail in Appendices |0 and lUl 
We neglect the k-dependence in ek,r occurring in the function Ft. Within this approximation, the functions Fi, Ft 
do not depend on the azimuthal angle (p', they only enter the k integral through the polar angle 9' . The 9' integral 
can be transformed into an integral over p. This p integral is performed using the formalism presented in Ref. |l9l | . 
To subleading order, we may set fc ~ q ~ /x, which allows us, for all phases we consider, to write the result of the 
azimuthal integral in the following power series of p, 



d^'f%^k,c{)^asU + 4^+vl^ ■ (81b) 
The coefficients as are chosen such that they add up to one, cf. Eq. In all cases we consider, we find 

a. = '"'^ , (82) 

and there are at most two nonzero coefficients oi, 02. In general, the coefficients and 77^'* depend on the polar angle 
6 of the vector q (cos6' — q • e^). However, in most of the cases we consider, they are constant. The only phase with 
angular-dependent coefficients rf'* is the mixed polar phase, while the only phase with angular-dependent coefficients 
Os is the transverse A phase. 

With ril = tiq = 2/3, which holds for all phases we consider, we arrive at (for details of the p integral, see Ref. |19j) 



{hhc^k^r= r^e^f d(g-M) V^.^^ln , , (83) 




with b as defined in Eq. ®. The value of d can be determined from the coefficients in Eqs. (|81|l . 

d = -^ H+V2 + 2{vi + vl)] ■ (84) 
Vo 

Consequently, the result of the k integral in the gap equation can be obtained by computing the coefficients of the 
power series in p, since the p integral is generic for all phases. This result is very similar to that of Ref. 0| , where 
the q integral on the right-hand side of the gap equation could be performed directly because of the trivial angular 
structure. In the present, more general, case, however, the q integral still has to be done. In Eq. (|83|l . this integral 
has been divided into its angular part and the integral over the modulus q, which can be restricted to an integral over 
a small region of size 26 around the Fermi surface, where 6 is much larger than the gap but much smaller than the 
chemical potential. This integral can now be done in the usual way, keeping the angular dependence in the excitation 
energy eq^s- We briefly repeat this calculation, which leads to Eq. (|91|l . more details for the case of two different gaps 
can be found in Ref. jl^ . With the approximation 

q,s - ek.r) In + B(ek,r - eq,s) In (85) 
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and the new variables 



Xr = g In 



k — ^ + ek,r / 

to subleading order the gap equation (|83|) transforms into 



= In 



2b fi e- 



(86) 



dys Us tanh 



Xo 



2T 



dys tanh 



2T 



where we defined 



g In 



2b 



/i e 



a^o = g In 



bfie-'^ 



(87) 



(88) 



with 00 = 0(2;*) — 0(2:^2) • In the case of an isotropic gap, 00 is the value of the gap at the Fermi surface. Note that 
Eq. H87|l corresponds to Eq. (79) of Ref. The difference to that equation, besides the angular integral, arises from 
the simplification we made above by omitting the wave function renormalization factor. 

While Eq. (|83|l in principle is a set of two equations, one for each excitation branch, labelled by the index r, we see 
from Eq. (|57|l . that, after the change of variables and neglecting sub-subleading contributions, one equation for the 
variable Xr is left. We can write this single equation for the renamed variable x in the form 



(4)k H^) ^ -r^l 




dy y tanh 



2T 



-ai X 



(89) 



The solution of this gap equation for T = is found in the usual way. Differentiating twice with respect to x yields 
a second-order differential equation. The two constants of the general solution are determined with the help of Eq. 
(|89|l and its first derivative at the point x ~ X2. One obtains 



ix) = 00 [cos(a;2 ~~ x) + ai (^2 — x\) sin(a;2 — x)] 



(90) 



An exchange of the indices 1 and 2 in this solution yields the same final result for 0o. In order to determine 0o, one 
inserts the solution H9(J|I into Eq. (|89|l and evaluates the equation at the point x = X2- Then, the integrals on the right- 
hand side of the equation are trivial. Using sm{a—j3) = sin a cos /3— cos a sin /3 and co&{a—f3) = cos a cos /3-|-sin a sin /3 
and the approximations sinxo — xq, cosxq ~ 1 (note that xo is parametrically of order 0{g)), we obtain 



= y —^£q[cosx2+ai{x2 — x*i)svtvx2] 



(91) 



where (ffc)j^ = {(-q)(i has been subtracted on both sides of the equation. With the definition for X2 in Eq. and the 
approximations 



sma;2 

we find (using oi -|- 02 = 1) 



g In 


( 




gin 


^ 2bii\ 


^00 j_ 



+ g (In \/\q.2 + d) sin 
- g (In + d) cos 




(92a) 
(92b) 




+ .9 (C + 1^) sin 




(93) 
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We use the abbreviations 



c 



d = 



(94) 



with C = ln(A°\A°^2)^^^j defined in Eq. iQ for isotropic gaps. With Eqs. (0 and (|H21l, we can write ( in the 
foUowing simple form, 



c 



1 (ni Ag^i In Ag,i + n2Ag,2 In Ag,2)q 

2 (?^iA„,i + n2A„,2)q 



(95) 



This expression shows that ( can be determined solely from the spectrum of the matrix L+, cf. Eqs. (jSHJi l|65(l . and 
(El- 

From Eq. (|93() we deduce that, to subleading order (reinserting the factor b'^), 



6o = 266Qe e fj, exp ( ~ ^ 



(96) 



This result, which is the main result of this paper, is the generalization of Eq. I^l. It says that, in the case of 
anisotropic gaps and/or an angular-dependent value of d, one has to replace the exponents d ^ d, C ^ C- It turns 
out that in most of the cases we consider, d is a constant number. Only in the mixed polar phase, see Sec. IIVI 
d ^ d. However, the modification of the other exponent, ^, plays an important and nontrivial role, especially for the 
determination of the ground state. 

In order to compare the result with the corresponding one in the theory of superfluid "^He, we consider the special 
case of only one gapped (but anisotropic) excitation branch, Xq 2 = 02 = 0. In this case. 



V (Aq.Oq 

and the quadratic mean of the gap is given by 

= 2 6&Q e^"* exp 



exp ( lnAg^i)q 



TT 



exp ( -2(-^9,ilnAg,i)q 



(97) 



(98) 



In this special case, the exponent involving the angular dependence of the gap, is exactly the same as in "^He, cf. Eq. 
(3.63) of Rcf. [l|. 



C. The critical temperature 



In order to determine the critical temperature Tc, we proceed similar to the above calculation of the gap: Starting 
from the gap equation (|5^ . we apply the method presented in Ref. before doing the drig-integral. The basic 
assumption is that the shape of the gap function does not change with temperature, i.e., we employ the following 
factorization of the temperature-dependent gap function. 



(x,r)~(/.(T) 



0(x,O) 



(99) 



where (f){T) = (j){x2,T) is the value of the gap at the Fermi surface at temperature T, and (/)(a;,0) is the zero- 
temperature gap function (j){x) computed in the last section. This ansatz, inserted into Eq. H89|) . yields at the Fermi 
surface 



(4 



17 dyytanh 



<y) 



0(2/, 0) 



dy tanh 



2T 



2T 



dy y tanh 



<y) 



2T 



0(2/, 0) 



(100) 
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where the second integral in Eq. (|89|l has been divided into two integrals: One running from xq to x^, with x^, = 
X2 — g ln(2K), K ^ 1, and one running from x^ to xi,- In Ref. [Tsj it is shown how the integrals in the angular brackets 
on the right-hand side of Eq. H 100(1 are evaluated. The result is 



0= l^e.in''"'^' 



47r V ttT, 



^JKaK^) ■ (101) 



This leads to 



^ = -e^ . (102) 

Note that, although the ratio Tc/ipo depends on the constant C, the absolute value for Tc does not, since the factor 
e~'' of the gap (f>o cancels the factor e'" on the right-hand side of Eq. p02|) . Consequently, in units of the critical 

temperature in the 2SC phase, T^^'^, the critical temperature only depends on the constant d, Tc/T^^'^ — e^'^. 

We recover as a special case the result for two isotropic gaps, Eq. IjlOl) . For the case of only one gapped excitation 
branch, and using Eq. (|97|l . the ratio between the critical temperature and the quadratic mean of the gap at the Fermi 
surface for T = reads 

J^l^ expQ(A,,i lnA,,i)q) . (103) 

Consequently, we find a nontrivial modification of the BCS relation — O.570o also in the case of a single, but 
anisotropic gap. Obviously, the BCS relation is recovered from Eq. H103|l for the case of a constant gap, since in this 
case Aq 1 = 1. Note that the result l)l()3|l is identical to the one in the theory of superfluid '^He, cf. Eq. (3.63) of Ref. 

D. The pressure 

The starting point for the calculation of the pressure is the effective action, Eq. ((45(1 . In the two-loop approximation, 
we can write 

r2[A,5] = iTr(S5) . (104) 



Then, making use of the Dyson- Schwinger equation 1(46 b() . the fermionic part of the effective potential at the stationary 
point can be written as 

T[S] - iTrln5-i - ^Tr (l - S^^S) . (105) 

In order to evaluate the first term on the right-hand side of this equation, we use the identity Trln S^^ — lndet5~^ 
and the fact that for arbitrary matrices A, B, C, and an invertible matrix V, 

det (^■^ ^det{AV- BV-^CV) . (106) 

Then, making use of — Sq^ + E and Eqs. 1(48(1 and ((50(1 . the trace over Nambu-Gor'kov space yields 

^Trln5-i = ^Trln {([G+]-i + S+) {[G^]-' + E") - {[G^]-' + S")-! $+ {[G^]-' + S")} . (107) 
In order to proceed, we set Z{ko) ~ 1 in the fermion propagator G^, cf. its definition in Eq. 1(61(1 . Using the identity 

mr'iG^r' = Y.l''l - in- ekf]At , (108) 



we find 

2 

1 
2 



^Trln5-i = ixr In ^ [fc^ - (^ - efc)^ - 02^+] a, 

e 

- J E E Tr[^^ .A-^] In [k^, - (el^f] . (109) 



e,r K 
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After performing the Matsubara sum, this expression reads 



2 T hi 



1 + cxp ( — ^ 



(110) 



where the trace now runs only over color, flavor, and Dirac space. 

The Nambu-Gor'kov trace for the second term on the right-hand side of Eq. (|105|l is easily performed with the 
definitions (|48() and (|51|l . Note that the anomalous propagators S^, occuring in the full quark propagator S, do not 
enter the result. One obtains 



iTr(i - s^'s) = -i E E nKAl + KA 



e,r K 

which becomes, after performing the Matsubara sum, 
^ j vn - ,s„ \st = — > / ■ 



tanh 



2T 



(111) 



(112) 



Using Eq. (fZ^ . the final result for the pressure p = —Vcs = yT, obtained by putting together Eqs. piU|l and (|112|1 . 



IS 



P = 



d^k 
(2^ 



1 + cxp ( — ^ 



2eS;„ 



tanh^ . (113) 



In the following, we restrict ourselves to the zero-temperature case, T = 0. Furthermore, we neglect the antiparticle 
gap and thus denote 4> = (j)^, as above. In this case, Eq. pi3|) becomes 



P 



d^k 



■ rir 



2e 



(114) 



k,r 



In order to evaluate the integral over the absolute value of the quark momentum, we assume the gap function to be 
constant in a small region around the Fermi surface of size 26, (pk^r = 4'o, and zero elsewhere. Moreover, we use the 
integrals 



Jo 



d^ 



1 



(115) 



and 



dk fc2(4o-rcj^.o 



e,.n - 2k) = -fi'^ , 



(116) 



where the vacuum energy 2k has been subtracted and e^.p = \k — e/i|. Neglecting terms of order 0((^g), we find 



P 



487r 



(117) 



where we denote the difference of the superconducting phase to that of the normal phase by 

,2 



Ap : 



16 TT^ 



E'^r^r- 



(118) 



As expected from physical intuition and from the theory of "^He, cf. Eq. the favored phase is determined by the 
largest condensation energy Ap, which is proportional to the sum of the angular averages of the squares of the gap, 
weighted with the corresponding degeneracies Ur of the excitation branches. (Note that, since 0^, = ^ (-^fc.r)^ 4'Q is 



—2 



the quadratic mean of the gap, 0^ = (Afe.r)j^ (f>o is the angular average of the square of the gap, not the square of the 
angular average, as the notation might suggest.) 

As a special case, one finds for the 2SC phase, counting Dirac, color, and flavor degrees of freedom. 



Ap2SC = 



2 7r2 



(119) 



This is in accordance with Ref. '30'| , when one identifies the pressure with the negative value of the effective potential 
at the global minimum. 
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Polar phase 




Mi, = Ji [a fca +/37±,3(k)] 




= Ji + (a' - cos" e] 




Afc.i =f3' + - IS"") cos' 9 (m = 8) , Afe.a = (na = 4) 








longitudinal 


mixed 


transverse 




0082 6* (8), (4) 


1/2 (8), (4) 


sin^e* (8), (4) 




1, 


1, 


1, 


d 


6 


5 


9/2 


c 


-1/3 


-In^ 


In 2 - 5/6 



TABLE L Relevant quantities for the polar phase. 



IV. DISCUSSION OF THE POLAR, PLANAR, A, AND CSL PHASES 

In this section, we use the general results of the previous sections for a discussion of the physical properties of 
certain phases in a spin-one color superconductor. In particular, we determine the ground state at zero temperature. 
We focus on the four "inert" phases presented in Sec.^ the polar, planar. A, and CSL phases, see Fig.^ 

The common structure of all spin-one phases is given by the matrix Mk, that determines the color and Dirac 
structure of the gap matrix $+(i4r), Eq. (|58|l . The most general form of this matrix has been introduced in Eq. H17|l . 
Specifying the spin-triplet structure Kj, we write the matrix as 



Ml 



akj -I- /37j.j(k) 



(120) 



The first term in angular brackets, proportional to the j-th component of the fermion momentum unit vector kj, 
describes pairing of quarks with the same chirality, since it commutes with the chirality projector Vr,e = (1 i 75)/2. 
The second one, proportional to 



7± J (k) = 7j - A;j7 • k , j = 1, 2, 3 , 



(121) 



corresponds to pairing of quarks of opposite chirality, since commuting this term with the chirality projector flips the 
sign of chirality. In the above ansatz for the gap matrix we allow for a general linear combination of these two terms, 
determined by the real coefficients a and (3 with 



-P^ = 1 



(122) 



In Refs. the special cases (a, /?) — (1, 0) and (a, /3) — (0, 1) were termed "longitudinal" and "transverse" gaps, 

respectively. We shall also use these terms in the following. (In Ref. fl2|, the LL and RR gaps correspond to the 
longitudinal and the LR and RL gaps to the transverse gaps.) The reason why both cases can be studied separately 
is that a purely longitudinal gap matrix on the right-hand side of the gap equation does not induce a transverse gap 
on the left-hand side and vice versa. More precisely, inserting the matrix A^k from Eq. I|12U|1 with f3 — into the 
anomalous propagator from Eq. H7U|I . and the result into the right-hand side of the gap equation H5()|) . we realize 
that the Dirac structure still commutes with 75 and thus preserves chirality. The analogous argument holds for the 
transverse gap, a = 0. For the case of an equal admixture of longitudinal and transverse gaps, i.e., a = /3 = l/V^, let 
us use the term "mixed" gap. For the physical results, we focus exclusively on the longitudinal, mixed, and transverse 
gaps. 

We summarize the relevant quantities for the different phases in Tables III - II VI The physically important results 
are collected in Tables M - IVlTl and Fig.EJ 

Let us first comment on Tables HI - II VI For the planar and A phases, we use the abbreviations 



^1/2 
B 
Z 



[a^ - /32) k^ , 



P^a fc2 7-L,i(k) - A:i7±,2(k) -P 7±,i(k) 7±,2(k) - ^1^2 | 



(123a) 
(123b) 
(123c) 
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Planar phase 




Mk = Ji[afci +/37±,i(k)] + J2[afc2 +/37±,2(k)] 




= J?^^ + J2A2 + {Ji,J2}B+ [Ji , J2] Z 




Xk,i=a''8m^e + f3\l + cos^e) (m = 8) , Afc,2 = (712 = 4) 




^^i = ikAfe.i , ^^2 = l-i^A.,i 




longitudinal 


mixed 


transverse 




sin^ei (8), (4) 


1 (8), (4) 


1 + cos2 6i (8), (4) 


Qr 


1, 


1, 


1, 


d 


6 


21/4 


9/2 


C 


In 2 - 5/6 





ln\/2 - 7/12 + 7r/8 



TABLE II: Relevant quantities for the planar phase. 





A phase 




Alk = J3{«(fci + ifc2) +/3[7±,i(k) +m,2(k)]} 




Li = Ji[{Ai+A2)±2i Z] 




Afc,i/2 = a^sin^6' + /3^(l+cos^6')±2/3 v'^^sin^fi' + ZJ^cos^e (jii = n2 = 4) , 

Afc,3 = (ns = 4) 








longitudinal 


mixed 


transverse 


Xk,r (rir) 


(8), (4) 


2 (4), (8) 


(1 + |cosei)2 (4), (1- |cose|)2 (4), (4) 




1, 


1, 


1/2 + 1 cos 611/(1 + cos^6l), l/2-|cos6l|/(l+cos^6l), 


d 


6 


21/4 


9/2 


c 


In 2 - 5/6 




In 2 - 1/3 



TABLE III: Relevant quantities for the A phase. 



The quantities A1/2, B, and Z are diagonal in color space. But while A1/2 and B are scalars, Z is a nontrivial 4x4 
matrix in Dirac space. One can verify the relation 

= B^ - A1A2 . (124) 

Moreover, in the table for the planar phase, we denote the anticommutator by {— , — }. The angle 6, used in Tables|l] 
- mil is the angle between the quark momentum k and the z-axis. The indices a,b < 3 in the second line of Table HVI 
are color indices. 

The specific form of the matrices A^k is obtained by inserting the respective order parameters A, cf. Fig.^ into 
the definition (|120|l . The matrices are given by the definition H63() . The calculation of their eigenvalues Xk,r and 
corresponding degeneracies Ur is presented in Appendix^ for the planar and A phases. The same method can be 





CSL phase 




A1k = J - [ak + /37x(k)] 




(L±)„b = {a^+2p^)Sab - [Qfc6 + /37±,i>(k)][aL -/37±,a(k)] 




Afc,i/2 = ^a' + 2(3'^ ± ia^a^ + 8/3^ (m = nz = 4) , Xk,3 = a' {m = 4) 




p± _tt3 ^i-^k.s 




longitudinal 


mixed 


transverse 


Xk.r (rir) 


1 (8), (4) 


2 (4), 1/2 (8) 


2 (8), (4) 


Qr 


1, 


2/3, 1/3 


1, 


d 


6 


5 


9/2 


c 





lnv/2173 


In a/2 



TABLE IV: Relevant quantities for the CSL phase. 



20 



rrr / i2SC 

\/ Ak,r(P0/4'0 


longitudinal 


mixed 


transverse 


polar 




e-5 




planar 


sine|ie5/«e-« 


e-21/4 


Vl+cos^e-ije'/^^-^/^e-^/^ 


A 


sin 6*1 ie^/^e"'' 


e-21/4 


(l±lcose|)iei/3e-»''2 


CSL 


e-6 


2(-l±3)/6 g-5 


e-9/2 



TABLE V: Gap functions Xk,r4>o in units of the 2SC gap. Fig. |5| illustrates this table schematically. 





longitudinal 


mixed 


transverse 


polar 


e-« 




e-9/2 


planar 


e-« 


e-21/4 


e-9/2 


A 


e-e 


e-21/4 


e-9/2 


CSL 






e-9/2 



TABLE VL The critical temperature Tc in units of the critical temperature for the 2SC phase. 



applied to the other phases, see for instance Refs. |0, l^^- After determining the eigenvalues, the corresponding 
projectors are obtained with the help of Eq. All these quantities, A^k, ^k,r, rir, and V^^, are computed for 

arbitrary linear combinations of longitudinal and transverse gaps and are given in the first four lines of the Tables ^ 
- II VI The last four lines of these tables are devoted to the special cases of pure longitudinal and transverse gaps and 
the mixed gap. First, we present the eigenvalues Xk,r with the corresponding degeneracies, immediately deduced from 
the general ones. Then, we compute the quantities a^, see Eqs. (|81|l . and d, see Eqs. II84|) and l|94l) . The calculation of 
these quantities is more or less straightforward, but may turn out to be lengthy. It involves performing the color and 
Dirac traces in the quantities T^g 0^^ l) ^^'^ (1*> Q) "^^ll as the angular integration over k, see Eq. (|80|) . For the 
most complicated cases, the transverse polar, planar, and A phases, we present details of the calculation in Appendix 
lUl Furthermore, in Appendix IdI it is shown that d = d = 6 is a universal result for all longitudinal phases. A special 
case is the mixed polar phase, which is the only case where d turns out to be angular-dependent, d = 3(3 + cos^ ^')/2. 
This result has already been obtained in Refs. [I^jllal. However, it has not been realized that c? = 5, and not d, enters 
the value of the gap ^o- Finally, the constant C is straightforwardly obtained by using Eqs. and l|94|) . 

Let us now turn to the physically important results. Tables IVl - 1 VIII and Fig. [3 In Table Ivl we present the (angular- 
dependent) gap functions for zero temperature at the Fermi surface, as they occur in the quasiparticle energies, Eq. 
(|67|l . They involve the square root of the eigenvalue A^.^ and the factors e^'^ and e~^, as shown in the solution of the 
gap equation, cf. Eq. I|96|l . The results are easily computed using the results of the previous tables. The magnitude 
of the gaps are reduced compared to the gaps in the spin-zero phases by factors of the order ~ 2.4 • lO^'^ through 
g-9/2 ^ y ^ ■ 10~2^ Consequently, assuming the spin-zero gaps to be of the order of lOMeV, the spin-one gaps are of 
the order of 10 - 100 keV. 

The angular structure of the gap functions, given in Table El is illustrated in Fig. El In this fi gure, the magnitude 
of the gap at the Fermi surface is shown as a function of the polar angle d. None of the functions depends on the 
azimuthal angle i.e., all figures are symmetric with respect to rotations around the z-axis. First note that the 
gaps in the CSL phase are isotropic. Nevertheless, there is a hidden anisotropy also in this phase, since the residual 
symmetry is not the group of rotations in real space but the group of joint rotations in color and real space. Thus, 
analogous to the B phase in ■^He one might call this phase "pseudoisotropic" . Other analogies to "^He can be found. 



Ap/Ap2sc 


longitudinal 


mixed 


transverse 


polar 


\e^l'A e-i2~0.65e-^2 




ie5/^e-9~0.88e-3 


planar 


~0.88e-^2 


e-21/2 


2e7/6-V4e-9^0.98e-s 


A 


ie^/^ e-i2^0.88e-i2 


1 e-21/2 


i 6^/3 e~9~ 0.65 


CSL 




3.2-4/3^-10 





TABLE VII: Zero-temperature pressure Ap (= condensation energy) in units of the 2SC pressure. The results are illustrated 
in Fig.m 
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particularly in the first column of the figure, representing the longitudinal gaps. All structures correspond to their 
analogues in ^He. This is plausible, because the respective gap matrices do not involve the nontrivial Dirac part 
7j^(k), and thus can, in this respect, be considered as the nonrelativistic limit. Note that the longitudinal polar phase 
has a nodal line at the equator of the Fermi sphere, while the longitudinal planar and A phases have nodal points at 
the north and south pole of the Fermi sphere. All longitudinal phases have one gapped and one ungapped excitation 
branch with degeneracies eight and four, respectively, cf. Appendix^ Note that the degeneracies add up to 12, since 
the matrix is a 12 x 12 matrix, involving antiparticle degrees of freedom. Nevertheless, the physical degeneracies 
of the gapped branches have to be reduced by a factor 2 compared to Fig. 12 since the antiparticle gaps are negligibly 
small. 

All phases shown in the second and third columns of the figure, have no analogues in '^He, because 7j^(k) gives rise 
to a nontrivial Dirac structure. The mixed gaps all are isotropic. The transverse gaps, however, exhibit nontrivial 
angular structures. The transverse polar phase has point nodes, similar to the longitudinal planar and A phases. The 
transverse planar phase has an anisotropic gap, however, the gap vanishes nowhere. The transverse A phase has one 
ungapped excitation and two different gapped ones, each with a nontrivial angular structure. One of these structures 
has no nodes but minima at the equator of the Fermi sphere, while the other has point nodes at the north and the 
south pole. 

long. mixed transv. 



polar 



planar 



A 



CSL 



/ 

magnitude of 
the gap 



FIG. 2: Schematic representation of the gap functions ^ Ar,fc (j>o, given in Table ^ and condensation energies, given in Table 
IVllI All gap functions are symmetric with respect to rotations around the z-axis. The numbers correspond to the degeneracies 
of the respective excitation branches. Only the gapped excitations are shown. The condensation energies are drawn to scale in 
each column separately. 

The critical temperatures, presented in Table IvTI are obtained from Eq. Hl()2|l . As remarked below that equation, 
the absolute value of Tc depends only on the factor e~'^. Therefore, the results of Table IVTI are easy to interpret. The 
constant d assumes the value d = 6 for all longitudinal gaps, which therefore all have the same transition temperature, 




condensation 
energy 



degeneracy of 
excitation brancli/ 
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being of the order of lOkeV (again applying (pQ^'^ — 10 MeV). Since the mixed phases have different values for d, their 
transition temperatures differ, ranging in the order 30keV - 40keV. The transverse phases have the largest transition 
temperature, Tc ~ 60 keV. 

The most important physical results of this paper can be found in Table IVTll and, schematically, in Fig. 13 since 
they answer the question of the preferred state in a spin-one color superconductor at zero temperature. Crucial to 
this question is the factor e^^ in the gap function, which is strongly affected by the nontrivial angular structures in 
the gap equation, which was not realized in Rcf. 28], where a simple approximation for the angular integral had been 
used. The first obvious result, already indicated by the results of Ref. JJJ and independent of the constant (, is that 
the pressure of all transverse phases is larger than that of the mixed phases, which, in turn, is larger than that of the 
longitudinal phases. Therefore, the most interesting and relevant phases are the transverse ones. Nevertheless, let us 
mention, that in the results for the longitudinal phases we recover Eq. Q, i.e., the ratio between the condensation 
energies of the longitudinal A and CSL phases is identical to the corresponding one of the A and B phases in '^He. For 
the transverse phases, the result cannot be anticipated from the theory of '^He. We find that the condensation energy 
of the phases with nodal lines or points of the gap, i.e., the polar and A phases, is smaller than that of the phases 
without nodes. However, note that the polar phase has a larger condensation energy than the A phase, although it 
has more excitation branches with nodes in the gap function (eight compared to four). The two phases without zeros 
in the gap function, i.e., the planar and CSL phases, differ only by two percent in their condensation energy. The 
isotropic transverse CSL phase has the largest condensation energy. 

Note that electric charge neutrality of the system does not affect the results, since we are considering a one-flavor 
system. This, of course, is in contrast to the systems where quarks of different flavors form Cooper pairs. In these 
systems, for instance in the (gapless) 2SC and CFL phases, charge neutrality plays an important role regarding the 
question of the ground state. The effect of the condition of color charge neutrality for a one- flavor system is less trivial. 
In the polar and A phases, Cooper pairs carry color charge anti-blue, i.e., only red and green quarks condense while 
the blue quarks remain unpaired. This pattern is known from the 2SC phase. It leads to a nonvanishing expectation 
value of the gluon field , which ensures color neutrality of the system. This is equivalent to introducing a color 

chemical potential /ig 6] which discriminates between the chemical potentials for paired and unpaired quarks. In the 
planar phase. Cooper pairs are formed by red/blue and green/blue quarks, i.e., they carry color charge anti-green or 
anti-red. Also in this case, one expects the color chemical potential for red and green quarks to differ from that for 
the blue quarks because of the residual color group SU{2). Consequently, quarks of different chemical potentials form 
Cooper pairs. In both cases, the condition of color neutrality might reduce the pressure. In contrast, the CSL phase 
automatically fulfills color neutrality, since there are equal numbers of Cooper pairs with colors anti-red, anti-green, 
and anti-blue. Therefore, the result that the transverse CSL phase has the largest pressure is insensitive with respect 
to electric and color neutrality conditions. 

V. SUMMARY AND OUTLOOK 

In this paper, we have investigated the possible color-superconducting phases in cold and dense quark matter, where 
quarks of the same fiavor form Cooper pairs, which leads to Cooper pairs with total spin one. It has been argued that 
these phases might be relevant in the interior of neutron stars, since a mismatch of Fermi surfaces, if too large, forbids 
pairing of quarks of different fiavors. Since the structure of the order parameter of the spin-one phases corresponds to 
a complex 3x3 matrix, a priori a multitude of phases seems to be possible. It has been the main goal of this paper to 
pick the order parameter that leads to the phase with the largest pressure, indicating that it is the favored spin-one 
color superconductor. In principle, this task is comparable to the three-fiavor case, where the CFL phase turns out 
to be the favored one. However, in that case, the argument that favors the CFL phase is more or less obvious and 
plausible, since it is the only possible phase in which all quasiparticles attain a gap in their excitation spectrum. The 
present paper shows that for the spin-one phases the argument is more subtle, since, in particular, it involves the 
angular structures of the gap functions. Angular-dependent gap functions are known from condensed-matter systems, 
for instance from superfiuid ^He. We have shown that a certain class of spin-one phases, named longitudinal phases 
(where quarks of the same chirality form Cooper pairs), reproduces some features of the phases in ^He, namely the 
angular structure and the relative condensation energies of the several phases. However, most of the investigated 
phases, and in particular the preferred transverse phases (where quarks of opposite chirality form Cooper pairs), have 
no analogues in ^He. 

The specific technical and physical results of the paper can be summarized as follows. In the systematic classification 
of order parameters, based on group-theoretical arguments, we have found four phases with uniquely defined order 
parameter matrices. These phases, the polar, planar. A, and CSL phases, all break color and rotational symmetries. 
From the theory of ^He we transfer the conclusion that these order parameters correspond to the inert phases, although 
the different group structure in principle demands for a more careful investigation of this statement. Nevertheless, we 
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have picked the above mentioned four phases for a more detailed discussion. Before evaluating these phases specifically, 
we have presented a general solution of the QCD gap equation. This treatment generalizes the one presented in Ref. 
[T^ . The main conclusion is that an angular dependence of the gap function gives rise to a nontrivial factor in 
the expression for the zero-temperature gap. The explicit result of the gap is 



00 = 2&6oe"'^e"S^ cxp ^-^j . (125) 



It generalizes Eq. ((SJ by the replacements 

The quantity dq is (half of) the sum of the eigenvalues \k,r of the matrix L^, cf. Eq. 1)78(1 . These eigenvalues contain 
the angular structure of the gaps. According to Eq. ((95|l . C, can be determined solely from the quasiparticle spectrum. 
The physical meaning of the replacement C C can be illustrated as follows. Consider an isotropic gap function, 0, 
and an anisotropic one, say | sin^j </>, where 9 is the angle between the quark momentum and one fixed spatial axis. 
One would immediately conclude that the condensation energy of the former one is larger, since the latter assumes 
its maximum value 4> only for 9 = 7r/2, while it is reduced for all other angles. However, we have found that the latter 
one actually is equipped with the above mentioned factor, i.e., it reads | sin 6*1 (p. If e^'» > 1, the question of the 
preferred phase is nontrivial and depends on the special value of C. The general result presented here reproduces the 
special cases of one or two isotropic gaps. 

Furthermore, we have computed a general expression for the critical temperature, 

^^-e<. (127) 

This equation shows that the BCS relation Td 4>q ~ 0.57 is not only violated in the cases of a two-gap structure and a 
gapless color superconductor. It is also violated in the case of a single gapped excitation branch if the corresponding 
gap is anisotropic. In this case, the ratio between the critical temperature and the quadratic mean <j) of the gap is 
given by Eq. H103I) . which is exactly the same result as for "^He. 

Finally, a general expression of the pressure at zero temperature has been derived starting from the effective action, 
which can be obtained from the QCD partition function using the CJT formalism. In particular, the pressure contains 
the condensation energy (density), which is the difference of the pressure in the superconducting phase compared to 
the normal-conducting phase, 

^P=^2T.^rt, (128) 

r 

Physically, this result is plausible. It is equivalent to counting the gapped branches and weight them with the angular 
average of the square of the corresponding gaps. Again, our general result reproduces the well-known results for the 
spin-zero gaps. 

These general results for the gap, the critical temperature, and the pressure have been applied to the above 
mentioned four spin-one phases. The genuinely new results concern all phases with nontrivial angular structure. 
These are the longitudinal and transverse polar, planar, and A phases, while the "pseudoisotropic" CSL phase already 
has been discussed in Refs. has been shown that several of the phases exhibit nodal points or lines of the 

gap, similar to '^He and several high- Tc superconductors in condensed-matter physics 33] . The magnitude of the gap 

is dominated by the factor e^"^ which varies from for longitudinal gaps to e~^/^ for transverse gaps. These factors 
have already been found in Refs. 0, 0|. The consequence is that the spin-one gaps are smaller by 2 - 3 orders of 
magnitude than the spin-zero gaps, which renders them of the order of 10 - 100 keV. These numbers are not changed 
essentially by the factor which, in all cases, is of order one. Nevertheless, this factor is decisive for the question 
of the favored phase, which can be found among the transverse phases. Since the pressure of all transverse phases 
contains the factor e^^ and since the number of the gapped excitation branches in all transverse phases equals 8, the 
preferred phase is determined by the specific factors originating from the angular structure. The result is that the 
transverse CSL phase has the largest pressure. This phase has an isotropic energy gap, contrary to all other transverse 
phases. It has been discussed that this result remains unchanged under the condition of overall color charge neutrality. 

Let us finally add some remarks regarding astrophysical consequences. As a first result, we note that the existence 
of a spin-one color superconductor in the interior of a neutron star is not ruled out by the temperatures, which, in the 
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case of old neutron stars, is in the range of keV. In order to find consequences of a spin-one color-superconducting core 
for the observables of the neutron star, it seems to be very promising to study the interplay of the superconductor 
with the magnetic field of the star. In Refs. 0,113 it has been discussed in detail that the mixed polar as well as the 
mixed CSL phase exhibit an electromagnetic Meissner effect, contrary to the spin-zero phases. From these results one 
concludes that it is very likely that also the other spin-one phases expel electromagnetic fields, although an explicit 
calculation has not yet been done. For the transverse CSL phase, this statement is obvious due to its symmetry 
breaking pattern. It might be interesting to investigate whether the angular structure, in particular the nodal points 
or lines, leads to anisotropics regarding external magnetic fields, for instance to angular-dependent penetration depths. 
Furthermore, external magnetic fields should be included into the discussion of the favored phase. As we know from 
•^He, above a certain threshold for the external magnetic field, the B phase is no longer the ground state of the system. 
Instead, the anisotropic A phase becomes favored. Therefore, it remains a physically important project for the future 
to investigate if the CSL phase is still preferred in the presence of a magnetic field. 

Another measurable property of a neutron star which is likely to be affected by a color-superconducting core is its 
temperature, or, more precisely, its cooling curve (= temperature as a function of time) Recently, it has been 

argued, that, unlike a spin-zero color superconductor, a spin-one color superconductor could explain the observed 
cooling curves 36]. The main ingredient of this conclusion has been the magnitude of the gap, i.e., the keV-gap 
compared to the MeV-gap of the 2SC or CFL phases. The physical mechanism behind this cooling curve is the 
emission of neutrinos, which dominates the cooling of the star after the first few seconds after its creation. Since 
one of the involved processes of neutrino emission is the Urea process, which requires the breaking of a Cooper pair, 
it might be interesting to take into account not only the effect of the magnitude but also of the nodal structure of 
the gap. Moreover, the nodes of the gap function would certainly affect the specific heat of the system. Note for 
instance that in the A phase of '^He (which has the same nodal structure as the longitudinal A phase in a spin-one 
color superconductor), the specific heat depends on temperature according to a power law, while in the B phase the 
temperature dependence of the specific heat shows an exponential behavior. Similar effects can be expected for a 
spin-one color superconductor and have to be included into a careful discussion of the cooling behavior of a neutron 
star. 
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APPENDIX A: ORDER PARAMETERS CORRESPONDING TO RESIDUAL GROUPS H = [/(I) x H' 

In this appendix, we evaluate the invariance equation (|25|l . This means that we determine all possible residual 
subgroups of the form H = U{1) x H' and the corresponding order parameters, using Eq. H24|) as an ansatz for the 
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generator for the residual U{1). The invariance condition can be written as a system of nine equations, 

as 



2c j All +i63Ai2 
2c] Ai2 - ifogAii 
2c] Ai3 



2^3 

2c] A21 +i&3A22 



2V3 
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2V3 
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+ 2c^ A31 +i63A32 = 
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(Alb) 
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(Aid) 
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The corresponding coefficient matrix A exhibits a bfock structure and the determinant thus factorizes into four 
sub-determinants. Therefore, we have to consider the equation 



= det A = det Ai detA2 det^s det A4 , 



where 



detAi 

det^a 

detA3 
detA4 



2%/3 
as , 



2c 



2V3 



2c 



as 
V3 



2c 



as 
V3 



2c . 



(A2) 

(A3a) 

(A3b) 

(A3c) 
(A3d) 



Now, one can systematicahy hst ah possibihties that yield a zero determinant of the coefficient matrix and thus allow 
for a nonzero order parameter. 

1. det^i = 0. 

Here we distinguish between the cases (i) where the two terms in the angular brackets of Eq. ljA3ap cancel each 
other and (ii) where they separately vanish. 

(i) as, c arbitrary, 63 = — as/(2V3) + 2c. 

Inserting these conditions for the coefficients into Eqs. HA1|I . one obtains for the order parameter matrix 



1 



Ai iAi 



A = - A2 iA2 




(A4) 



where the factor 1/N with N = (2|Aip + 2|A2p)^^^ accounts for the normalization (|26|l . In this case, the 
order parameter contains two independent parameters Ai and A2. From this form of the order parameter, 
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we can now determine the group H' in Eq. (|22() . Inserting A into Eq. (|21(l and using the fact that the 
parameters Ai, A2 are independent of each other, one obtains the conditions 



fli = . . . = 07 = 61 = 62 = 



2V3 



as 



63 - 2c = . 



Consequently, 



H = U[l) X C/(l) 



(A5) 



(A6) 



since a vanishing coefficient in Eq. (|A5|) translates to a "broken dimension" of G. For instance, ai = means 
that Ti does not occur in the generators of H , etc. The dimensions of the residual Lie group can be counted 
with the help of the number of the conditions for the coefficients. Since dim G = dim Gi +dim G2 +dim G3 = 
8+3+1 — 12, and the number of conditions in Eqs. ljA5|l is 10, we conclude dim H = 2, which is in agreement 
with Eq. (|A6I) . Or, in other words, there is an additional U{1), i.e., H' = U{1) in Eq. H22() because the 
equation relating the three coefhcients as, 63, c allows for two linearly independent generators U and V 
which are linear combinations of the generators Tg, J3, 1. Note that U and V are not uniquely determined. 
One possible choice is 



C/-Ts- 



1 



2V3 



J3 



y = 2J3 + 1 



(A7) 



Different order parameters are obtained from two subcases: 

First, one can impose the additional relation c = —as/{2^/3) between the two coefficients that have been 
arbitrary above. Then, ^3 — 



These two conditions yield 
/ Ai iAi 



A = 



1 

N 



A2 iA2 
V A3 



(A8) 



where TV 
a 



(2|Ai| 



m, bn, c, which leaves a subgroup 



A3p)^/^. In this case, Eq. ()21fl leads to 11 conditions for the coefficients 

(A9) 



H = [/(I) , 

generated by a linear combination of generators of all three original subgroups Gi, G2, G3, 

1 



U ■ 



2V5 



(AlO) 



Second, one can set one of the coefficients Og, c to zero. The condition c — does not yield a new case. 
But ag — 0, and consequently 63 = 2c, has to be treated separately. In this case, Eqs. IjAljl yield 



A = l 

N 



where N = (2|AiP + 2IA2I 



generated by 



(ii) c = ag/(4y3), 53 = 0. 
Here, one obtains 



^ Ai iAi 
A2 iA2 
VA3 iAs 

2|A3p)^/^. The residual group is given by 
H^U{1) , 

[/ = 2J3 + 1 . 



(All) 



(A12) 



(A13) 



1 

N 



Ai A2 A3 
A4 A5 Ae 




(A14) 



where iV = eti|A.|2)i/2. A gain, the residual group is one-dimensional, 



generated by 



u = n 



2. det^2 = 0. 

This determinant vanishes in the following cases: 

(i) &3 arbitrary, c — as/{4:\/3). 
With Eqs. (jAip . one obtains 



/ Ai 
Aa 
^0 



where N = ([Aip + jAap)^/^. Inserting A into Eq. ^ yields 



ai 



a? = &i = &2 = , c : 



4^3 



aa 



As for the order parameter l|A4|l . the residual group is given by 

H ^ U{1) X C/(l) . 
However, the corresponding generators differ from those in Eqs. (|A7|I . 



4V3 ' 



y - J: 



3 • 



(ii) &3 arbitrary, ag = c = 0. 
In this case, 



1 



Ai 



A=- A, , , 
A3 

where N = (jAip + jAjp + |A3|2)i/2, xhe residual group is 

H = Uil) , 

which is a subgroup of the spin group G2 = SU{2)j, since it is generated by 

C/ = J3 • 

3. det^3 = 0. 

(i) ag, c arbitrary, 63 — aa/\/3 + 2c. 
In this case, we find with Eqs. (|Af p . 



A„| 000 







1 i 
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This matrix differs from all previously discussed order parameters in that it is uniquely determined. It 
corresponds to the A phase. Here, as in all cases above, we omitted a possible phase factor which could 
multiply A without violating the normalization. Inserting A into Eqs. H21(l yields the following relations, 



a4 = . . . = 07 = 61 = 62 = , ^pas — 63 + 2c = 

v3 



(A25) 



Consequently, besides the relation between og, 63, and c, there are only 6 additional conditions. Thus, the 
dimension of the residual group is 12—7 = 5. We obtain 



H = SU{2) X U{1) X C/(l) , 



(A26) 



where SU{2) is generated by Ti, T2, T3, and thus is a subgroup of the color gauge group Gi = SU{3)c 
For the generators of the two [/(l)'s one can choose 



2V3 



(A27) 



As in case l.(i), there is a subcase that produces an additional order parameter. Namely, if we require the 
condition c = ag/(4\/3), which yields 63 — %/3a8/2, we obtain 



/ A2 
A3 
y Ai lAi 



(A28) 



where N = (2|Ai|2 + jAap 
Consequently, 



with the generator 



|A3|2)i/2, From Eqs. we conclude that all other coefRcients vanish. 



(ii) &3 = 0, c = -a8/(2V3). 
With Eqs. (HU one obtains 

^ / 
A=- 
\Ai A2 A3 

where N is defined as in Eq. l|MT1l . From Eq. ^ we conclude in this case 



(A29) 



(A30) 



(A31) 



04 



aj = bi 







Hence, the residual group is 



generated by Ti, T2, T^, and 



2V3 



H = SU{2) X U{1) 



(A32) 



(ASS) 



4. det^4 = 0. 

There are two cases in which det^4 = 0: 



2x/3 



(AS4) 
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(i) &3 arbitrary, c = —ag,/{2^/?>). 
These relations lead to 



(A35) 

As in the A phase, Eq. ljA24p . the order parameter is uniquely determined. It describes the polar phase. 
Inserting the order parameter of the polar phase into the invariance condition, Eq. I)21|l. yields 

1 

a4 = . . . = 07 = 61 = 62 = , c = j=aB ■ (A36) 

2v 3 

Therefore, the residual group is 

H = SU{2) X [/(I) X U{1) (A37) 

with the generators Ti, T2, T3, and 

C/ = T8-^1, V = .h. (A38) 



Thus, the symmetry breaking pattern in the polar phase is similar to that in the A phase, cf. ljA27p . But 
while in the A phase both residual t/(l)'s are combinations of the original symmetries, in the polar phase, 
one of them is a subgroup of G2 = SU{2)j. 

(ii) 63 arbitrary, as = c — 0. 

This case is identical to case 2.(ii). 

APPENDIX B: EIGENVALUES OF L± IN THE PLANAR AND A PHASES 

1. Planar phase 

In order to compute the eigenvalues of the matrix = L^, given in Table (nj we use Eq. H124|) and the 
(anti)commutation properties of the color matrices {Ji)jk ~ ^i^ijk to obtain {L^Y = (^1 +^2)ik- Therefore, 

{L^r^iA.+A^T-'Lt . (Bl) 
This simple relation for the n-th power of can be used to find the roots of the equation 



det(A - lJ) = , (B2) 



which yield the eigenvalues A of L J . Since 



det(A-Lj) EE exp[Tr In (A- , (B3) 



we consider the trace of the logarithm of A — , 



Trln(A-L^) = In A Tr 1 + Tr In 1 - ^ = In A Tr 1 - ^ - A"" Tr . (B4) 



A / ^-^ n 

^ n—l 

The second term on the right-hand side can be evaluated using Eq. With TrL^ =8(^1+^2) one obtains 

det(A-L±) = A^[A-(Ai+A2)]^ (B5) 

which yields the eigenvalues Ai + A-i and with degeneracies 8 and 4, respectively. With the definition H123a|l . this 
confirms the entries for Afc r in Table HII 
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2. A phase 

In order to derive an expression for the n-th power of the matrix , given in Table IIIII we note that 

[L^f = 2(Ai + A2)L^ - [(Ai - ^2)2 + 4^2] J| , (B6) 

where Eq. 11241) has been used. Since A1/2 and B are scalars and J^L^ — L^, we have 

= a„i± + 6„ J| , (B7) 

with real coefheients a„, 6„. Applying Eq. (|B6|I . one derives the following recursion relations for these coefficients, 

an+i = 2{Ai + A2) an + bn , bn+i = -[{Ai~A2f +4B^]an . (B8) 

With the ansatz of a power series, a„ — p", one obtains a quadratic equation for p, which has the two solutions 

Pi/2=Ai+A2±2^AiA2-B^ . (B9) 

Therefore, a„ is a linear combination of the n-th powers of these solutions, a„ — 771J3" + fy2P2 • The coefficients 771, 772 
can be determined from oi = 1, 02 = 2(Ai +742). One finds 

""~4V AiA2-i3^ ^^^ 4 VMr~B^ ^ ' 

With TrLj =8(^1+^2) this yields 

det(A - L±) = A4(A - pi)4(A - p2)^ . (Bll) 

Consequently, both and iThave the eigenvalues 0, pi, and p2, each with degeneracy 4, which confirms the 
corresponding entries in Table 



Iml 



Finally, we consider the projectors V^^, corresponding to the eigenvalues A. With 



7?± ^ -^k(-^k ^2/1) p± _1_T?± _'p± (-Q12] 

' k.1/2 ~ \ ^\ \ ^ ' ' k.3 ~ ' k.l ' k.2 ' l^^^J 



k,l/2 \ , /'X , \ , A ' k,3 ' k,l ' k,2 ' 

one obtains 



p± = - J2 1 ± z , (B13a) 



^k,3 - 1 - -^1 ■ (B13c) 
Since TrZ = Tr(Aj;.Z) = 0, Eq. ^ is obviously fulfilled also in the A phase. 

APPENDIX C: TRACES AND ANGULAR INTEGRALS FOR THE TRANSVERSE PHASES 

In this appendix, we present the technical details of the solution of the gap equation for the transverse polar, 
planar, and A phases. The transverse CSL phase is trivial with respect to the angular structure and has already been 
discussed in Ref. . We compute the color and Dirac traces and perform the k-integration of Eq. H80(l in order to 
determine the coefficients and 77^'* introduced in Eqs. (|81|l . These coefficients yield the constant d, cf. Eq. (|84|) . 

Let us define the following traces over Dirac space. 



^ Tr 
7e™" ^ Tr 



7M7±,™(q) Aq 7^7j_,„(k) A+J , (Cla) 
7M7-L,m(q) 70757 •qAq7i^7J-,n(k)A+ . (Clb) 
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We shall make use of the following results, 



00 



■mn 
ij 



'^00 '^00 



'5r?xri(l + Q ■ k) + kmkn + (jmQn Qmkn q ' k + qnkm j (C2a) 
2 -PmPn + q-nkm - Q ' k ((5mn - PmPn) + Pm{kn - (7„) P ' k 

-Pn(fcm - 9m) P • q+ ((^mn - 9m fcn) P • qP • k , (C2b) 

*(g3 + fc3)(l + q-k) , (C2c) 

-2z {[l-(p.k)2]g3 + [l-(p.q)2]fc3_(l_q.k)(p.k + p.q)p3} . (C2d) 



1. Transverse polar phase 

Using the matrix A^k from TableQlwith (a,/3) = (0, 1), we find 



Tr 



4(1-^3') , 



and thus 



^oo(k,q) = 



1 + q • k - ^3^3 (1 - q • k) - ^2 _ ^1 



where the color trace Tr[T^ JsTaJs] = —4/3 has been used. Analogously, we find 
2 



^t'(k,q) = 



3(l-fc|) 



-93^3 + q - k + P3 (1 - q - k) - p3 (93 - ^3) (p • q - p • k) - (1 - 93^3) p • qp • k 



(C3) 



(C4) 



(C5) 



In order to obtain these results, Eqs. (|C2ap and ljC2bp with m = n = 3 have been employed. As expected, the 
ungapped excitation branch does not contribute to the gap equation, 7^Q(k, q) — 7^2(k, q) — 0, and thus ai — 1, 
a2 = 0. 

Next, the angular integral dflk has to be performed. To this end, we use the following frame (cf. left diagram in 
Fig. El : The order parameter in the polar phase picks a special direction in real space. By convention, we choose 
this direction to be parallel to the z-axis. Now there is another fixed direction, q, which we can choose to be in the 
ccz- plane, q = (sin0, 0, cos6'), where 9 is the polar angle. In order to perform the dflk integral, we use a frame with 
2;'-axis parallel to q. The original frame can be transformed into the new one with a rotation R{9) around the y-axis, 



m 



cos 9 0— sin t 

1 
sin 9 cos 9 



(C6) 



In the new frame, we have q' = (0,0, 1) and k' = (sin 6*' cos (^', sin 6*' sin (p', cos 0'). Therefore, before performing the 
dQk integral, we write k as 



k== i?-i(6') k' 



/ cos 9 sin 9' cos (p' + sin 9 cos 9' 

sin 9' sin tp' 
y — sin 9 sin 9' cos (p' + cos 9 cos ( 



(C7) 



This new frame is particularly convenient, since the new polar angle 9' is the angle between q and k, cos 0' = q • k. 
Therefore, the integration over 9' can be transformed into an integral over p via cos^' = (/c^ + —p'^)/{2kq), where 
p is the modulus of p = k — q. Note that the function on the right-hand side of Eq. H80() does not depend on the 
azimuthal angle ip' but only on p, and the function Ft only depends on p and 9, but not on ip' either (neglecting the 
k dependence in ek,r)- Consequently, we have to multiply T^q and T^} , given in Eqs. ljC4|l and ljC5|l . respectively, with 
the factor £k/£q = (1 — kl)/{l — write the results in terms of the new coordinates and integrate over (p' . After 
setting k ~ q ~ n, which is permissible to subleading order, the result yields the coefficients 77 in Eq. (|81|l . Then, the 
result of the p integral obtained by Eq. (|84|) , following Ref . [l^ . 
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q = z 





polar phase 



X q = z 

planar phase 



FIG. 3: Rotation of the coordinate system {x,y,z) {x',y',z') for the k-integration in the cases of the transverse polar and 
planar phases. 



The explicit results of this procedure are as follows. With 



2Tr 



^ / d(p' ki ^ q, q ■ k 



^ J^'^ d^' kf = qf (q • k)2 + i (1 - qf) [l - (q • k)^' 



we obtain 



where the approximation k q ^ has been implemented via the replacements q- k — > 1 — p^/(2/i^), p • k : 
p/(2/i). Now, we immediately compute d = 9/2, as listed in TablcQ] 



(C8a) 
(C8b) 

(C9) 

p 



2. Transverse planar phase 



In this case, we find from Table Ull 



Tr 



4(l + fc2) 



Furthermore, 



K.i = T—^ [-^I'^l - + ^2(1 - il) - {Ji, J2} 9192 + ^393 70757 ' q] 
+ 93 



where we have applied the identity 



7-L,i(q)7±,2(q) - qiq2 = iqaioi^i ■ q ■ 



(CIO) 



(Cll) 



(C12) 



With the help of Eqs. (jC10|l . HC11|I and the definitions ifT^ . fZ^ . we compute the quantities 7j3Q(k, q) and 7^^(k, q). 
As in the polar phase, we find 7^Q(k, q) = 7^^(k, q) = 0, hence ai = 1, a2 = 0. For the gapped branch, s = 1, one 
obtains after taking the trace over color space, 



T^Uk,q) = - 



{ql g| - 2) Qll + (gj - g| - 2) Qf^ + q,q, {Qll + Q^i) + ^qs {TZH ~ Tlf^) 

6(l + fc|)(l + gi) 



(C13) 



33 



where the tensors Q and TZ are defined in Eqs. (jCl|) . In order to perform the traces over Dirac space, one makes use 
of the identities HC2() . We do not present the explicit results for Tqq (k, q) and 7^^ (k, q) since they are too lengthy. 
In order to perform the angular integration, we proceed in a similar way as discussed above for the transverse polar 
phase. The difference is that in the planar phase the order parameter does not point into a special direction, but is 
located in the xy-plane. Without loss of generality, we can assume q to be also in the xy-plane, q — (cos sin(p, 0) 
(cf. right diagram in Fig. 13). Note that this choice, in particular ~ 0, immediately shows that the term proportional 
to 75 on the right-hand side of Eq. ljCll|) yields no contribution to the integral. Again, we rotate the frame such that 
the z'-axis is parallel to q. This can be done with two successive rotations i?i and i?2(</5): First, we rotate the original 
frame by 7r/2 around the y-axis. 



Ri = 



0-1 



1 
^10 



(CM) 



Then, we rotate by tp around the new x'-axis. 



/ 1 

cos ip — sin if 
y sin if cos if 



(C15) 



As above, in the new frame, q' — (0,0, 1) and k' = (sin 0' cos </?', sin^' sin</j', cos6''), and hence 



— sin 9' sin ip' sin ip + cos 9' cos </? 
sin 6' sin tp' cos tp + cos 9' sin p 
— sin 9' cos (p' 



(C16) 



With ik/iq = (1 + ^3)7(1 + 93) and the approximation /c ~ g ~ ^, we obtain 



From these results we conclude d = 9/2. 



2 lp2 I pi 

3 " 3^ ^ 24^ 



(C17) 



3. Transverse A phase 

In the transverse A phase, we have 



Tr 



Ail + kl), (C18) 



The projectors corresponding to the gapped excitations, ^qs- (|B13p . are 



q4/2 



where Eq. ljC12ll has been used. Inserting this projector and the respective quantities from Table UTTl into Eq. H72|l . 
one obtains after taking the trace over color space. 



Tr 



[7±,i(q) + m,2(q)] (1 ± |fT7o757 ' q) Aq 7.^ [7±,i(k) - m.2(k)] A+ 



T^/^ k,q = ^ , L , C20 

^ ' 6(l + fc|) ^ ' 

and 7^^j,(k, q) = 0. Again, the results for the quantities 7^o^(k, q) and 7J^'^(k, q), corresponding to the gapped 

branches, are complicated functions of k and q, and we do not present their explicit forms. 

As in the transverse planar phase, the order parameter for the transverse A phase points into a direction perpen- 
dicular to the z-axis. However, in order to perform the angular integral over 1^9', we cannot make use of the same 
choice of the frame, because the assumption of q being in the xy-plane, i.e., 93 = 0, would lead to a division by zero 
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according to Eq. (|C20|I . Instead, we allow for the most general form, q = (sin 6* cos (^s, sin sinews, cos 6'). In this case, 
the rotation of the original frame into the one with z'-axis parallel to q is given by 



cos t) cos (fi cos sm ip — sm ( 
— siniy9 cosiy9 I • (^21) 



Rie,ip) = 

y sin cos if sin 9 sin ip cos 
Consequently, the dflk integral has to be performed with 

cos Lf cos 9 sin 6' cos cp' — sin ip sin 9' sin p' + cos p sin 9 cos 9' 
k = i?^ ^ (6', (ys) k' = I sin p cos 9 sin 9' cos p' + cos p sin 6*' sin p' + sin sin 9 cos | , (C22) 

— sin 9 sin 9' cos i^' + cos 9 cos ^' 

where k' = (sin0' cos(y3', sin6'' sin(y9', cosf?') has been employed. 
With k q ^ one obtains 



27r Jo °" ^ '^^ 27r Jo 2 ^3 3 /i^ 24 /iV ^ 1 + cos2 e 



The term ± 2 | cos0| / (l+cos'^ 9) on the right-hand side of this equation gives rise to a difference between the quantities 
TqIj j(k, q) and T^q j(k, q). The origin of this difference is the term proportional to 75 on the right-hand side of Eq. 
(jC19ll . As a consequence, one cannot choose the coefficients as such that they are constants, cf. Eqs. (|81|l . The 
transverse A phase is the only phase we consider, in which these coefficients depend on 0, the angle between q and 
the z-axis. We find 

1 , Icos^l 1 Icos^l 

"^ = 2 + Tlw^' "^ = 2-TT^' (^''^ 

However, the constant d is identical to all other transverse phases, d — 9/2. 



APPENDIX D: PROVING d = 6 FOR ARBITRARY LONGITUDINAL GAPS 

In this appendix, we prove d = 6 for any order parameter {— for any 3x3 matrix) A in the case of a longitudinal 
gap. The longitudinal case is particularly simple since the Dirac structure of the matrix A^k is trivial. 
Let us start with the matrix 

Mk = t^k-J, (Dl) 

where = (^k,!; v\^.2, i'k.3) is a 3- vector with 

3 

^;k,^ = ^AylJ, 1 = 1,2,3. (D2) 

i=i 

Then, 

{L+),j = vl % - Wkj Wk,* , (D3) 

where = • Vk. Now, with {L'^Y = Wk-^k '^'"-^k " ^''^k the eigenvalues of are easily found making use 
of the method shown in Appendix IBl One obtains 

Afc,i = vl (8-fold) , \ka - (4-fold) . (D4) 

Consequently, the fact that there is one gapped branch with degeneracy 8 and one ungapped branch with degeneracy 
4 is completely general, i.e., it is true for any order parameter in the longitudinal case. The corresponding projectors 
are given by 

= S., - , iVi,h - . (D5) 
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This leads to 

1 Vr, ■ Vt . _i „ , 2 Vr, ■ V 

3 



T,l{k,q)^-^^il + ci-k) , r/(k,q) = -^J-^(l-p.qp.k) , (D6) 



and T^Q^k, q) = 7^^(k, q) = 0, hence ai = 1, 02 = 0. The angular integration is done as explained for the transverse 
A phase in the previous appendix, i.e., the rotation given by Eq. (|C21I) and the respective expressions for the vectors 
k and q are used. 

With = Afc^i/Ag,! = v^./v'^ we obtain 



where Eq. IjCSap has been used. Therefore, we find with fc ~ g ~ /i, 

2nJo ^ i, 2ttJo ^Iq*"-'^' 3 2 12 ^i* ' ^ ' 

which, making use of the definitions H81|l and (|84|) . proves the universal result d — 6 for longitudinal gaps. Since d — d 
when d is constant, cf. definition H94() . we conclude d~Q. 
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